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Abstract 

It is known that each symmetric stable distribution in is related 
to a norm on that makes M'^ embeddable in Lp{[0, 1]). In case of a 
multivariate Cauchy distribution the unit ball in this norm is the polar 
set to a convex set in M.'^ called a zonoid. This work interprets general 
stable laws using convex or star-shaped sets and exploits recent ad- 
vances in convex geometry in order to come up with new probabilistic 
results for multivariate stable distributions. In particular, it provides 
expressions for moments of the Euclidean norm of a stable vector, 
mixed moments and various integrals of the density function. It is 
shown how to use geometric inequalities in order to bound important 
parameters of stable laws. Furthermore, covariation, regression and 
orthogonality concepts for stable laws acquire geometric interpreta- 
tions. A similar collection of results is presented for one-sided stable 
laws. 
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1 Introduction 



Since P. Levy it is well known that the characteristic function of a symmetric 
stable law in R"^ can be represented as the exponential of the norm as (p{u) = 
e~ll"ll'', where p is the characteristic exponent of the stable law. It is also 
well known [20] that all norms which might appear in this representation 
make (M*^, || ■ ||) isometrically embeddable in the space Lp{[0, 1]) if p > 1. The 
corresponding result holds also for all p G (0, 2], see |26l Lemma 6.4]. 

Each norm in M.'^ gives rise to the corresponding unit ball F. In this 
paper we neglect the convexity property of the norm and use the term norm, 
where other works sometimes use the term gauge function. If the norm is 
convex, it can be realised as the support function = h{K,u) for the set 
K polar to the unit ball F, see Section |2l It is known from convex geometry 
(see IH] for the standard reference) that unit balls in spaces embeddable 
in Li([0, 1]) are exactly polar sets to zonoids. Recall that zonoids appear 
as Hausdorff limits of zonotopes, i.e. finite Minkowski sums of segments. 
In application to stable laws with characteristic exponent p = 1, we have 
that f{u) = e~^^^''^^ is a characteristic function (necessarily corresponding 
to the multivariate Cauchy distribution) if and only if is a zonoid. It is 
known that all planar centrally symmetric convex sets are zonoids, while this 
is no longer the case in dimensions 3 and more. This corresponds to a result 
of Ferguson [S], who showed that the dependency structure of symmetric 
bivariate Cauchy distributions can be described using any norm in M?, while 
such representation is no longer possible for all norms in dimensions three 
and more. 

It is explained in Section [3] that the correspondence between norms and 
stable laws can be extended to include all symmetric stable laws by repre- 
senting their characteristic functions using (possibly non-convex) norms. The 
unit ball F in the corresponding norm is a star-shaped set called the star body 
associated with the stable law. Thus, each symmetric stable distribution is 
uniquely determined by a star body F and the value of the characteristic 
exponent a and the paper aims to show that this geometric interpretation 
is useful in the studies of multivariate stable laws. In particular, F is an 
ellipsoid if and only if the underlying distribution is sub-Gaussian. Section H] 
shows that if the characteristic exponent p is at least one, then the norms 
are convex and the support function representation is also possible using the 
associated zonoid K being the polar set to F. This associated zonoid K 
is called Lp-zonoid, since K can be represented as the Hausdorff limit of a 
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power sum of segments, where the support function of K is the p-mean of 
the support functions of the summands. 

While the general correspondence between zonoids and stable laws is 
well understood, this paper concentrates on further relationships between 
probabilistic aspects of stable laws and geometric properties of associated 
star-shaped and convex sets. The core of the paper begins in Section where 
it is shown how to relate the value of the probability density function / of the 
symmetric stable law at zero to the volume of the associated star body F. It 
is shown how derivatives of / at the origin are related to further geometric 
properties of F, in particular to certain ellipsoids associated with F. It also 
provides an expression for the Renyi entropy of symmetric stable laws. Using 
geometric results on approximation of convex sets with ellipsoids. Section O 
ends up with a result that gives an estimate for the quality of approximation 
of a symmetric stable law with a sub-Gaussian one. 

Section [6] uses the Fourier analysis for generalised functions together with 
the geometric representation of the characteristic function in order to com- 
pute a number of important characteristics of symmetric stable laws. These 
characteristics include the moments of the norm of a stable random vec- 
tor, which previously were known only in the isotropic case, mixed moments 
of (possibly signed) powers of the coordinates, integrals of the density over 
subspaces, etc. This section also presents a number of inequalities for the 
moments and settles the equality cases. Finally, it clarifies a relationship 
between zonoids of stable laws and zonoids of random vectors studied in 

Section [7] deals with one-sided strictly stable laws supported by M!^. It 
first establishes a geometric characterisation of stable laws for power sums, 
which fill the gap between the arithmetic addition and the coordinatewise 
maximum scheme for random vectors. Note that relationships between max- 
stable random vectors and convex sets have been explored in [13] . It is known 
that max-stable distributions with unit Frechet marginals are exactly those 
having the cumulative distribution function F(u~^) = where is 

the vector composed of the reciprocals of u and is a max-zonoid, i.e. the 
expectation of a random crosspolytope. It is shown in Section [7] that stable 
laws for power sums (and also one-sided strictly stable laws) correspond to 
a new family of convex sets called Li(p)-zonoids. These sets appear as (set- 
valued) expectations of randomly rescaled £g-balls in W^, where q is reciprocal 
to p. Finally, it provides expressions for some moments of one-sided strictly 
stable laws. 
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The star bodies and zonoids associated with stable laws are determined 
by the spectral measures of stable laws. Section [S] shows that under quite 
general conditions, the spectral measures themselves admit a geometric in- 
terpretation as, e.g., surface area measures of further auxiliary convex sets 
called spectral bodies of stable laws. It shows how geometric inequalities can 
be used to relate volumes of the corresponding convex sets, and thereupon 
derive bounds for densities and moments of stable laws. 

The geometric interpretation of the covariation is given in Section [91 
It also discusses the regression problem for symmetric stable laws, in par- 
ticularly, the linearity property of multiple regression, which goes back to 
W. Blaschke's characterisation theorem for ellipsoids. 

Section dU] describes several operations with associated star bodies and 
zonoids and their probabilistic meaning. Using recent approximation results 
from convex geometry, it is proved that each symmetric stable law can be 
obtained as the limit for sums of sub-Gaussian laws. It also discusses optimi- 
sation ideas, which appear, e.g. in optimising a portfolio whose components 
have jointly stable distribution. 

Section [TT] discusses the concept of James orthogonality for symmetric 
stable random variables, which is also extended to define orthogonality of 
symmetric stable random vectors. 

This work attempts to highlight novel relationships between convex ge- 
ometry and the theory of stable distributions. Further developments are 
surely possible by invoking other recent results on isotropic bodies, geomet- 
ric concentration inequalities, or properties of convex sets in spaces of high 
(but finite) dimension. It is possible to apply the results to deal with finite- 
dimensional distributions of stable processes or work directly in a general 
infinite dimensional setting. 

2 Star bodies and convex sets 

A set F in is star-shaped if [0, C -F for each u E F. A closed bounded 
set F is called a star body if for every u E F the interval [0, u) is contained 
in the interior of F and the Minkowski functional (or the gauge function) of 
F defined by 

\\u\\f = inf{s > : M G sF} 
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is a continuous function of m G M . The set F can be recovered from its 
Minkowski functional by 

F = {u: \\u\\f < 1} , 

while the radial function 

Pf{u) = \\u\\p^ 

provides the polar coordinate representation of the boundary of F for u from 
the unit Euclidean sphere In the following we usually consider origin- 

symmetric star-shaped sets and call them centred in this case. If the star 
body F is centred and convex, then \\u\\f becomes a convex norm on W^. 
The ip-ball in M'^ is defined by 

B^ = {xeR'^: \\x\\p < 1}, 

where ||x||p = + ■ ■ --I- 1x^1^)"^/^ for p 7^ 0, i.e. ||x||p = ||a;||i7' with F = B^. 

If p G M \ {0}, the p-star sum of two star bodies Fi and F2 is defined by 
its Minkowski functional as 

\\^Kr+,F, = (Mf. + MfJ'^' ^ (2-1) 

This definition goes back to Firey |9j and was later investigated by Lutwak 
|33j . For p = — 1 we obtain the radial sum, i.e. the radial function of the 
result is the sum of two radial functions of the summands. Extended by the 
limit for p = —00, the p-sum yields the union Fi U F2 and the intersection 
Fi n F2 for p = 00. Note that (12.11) means that the Minkowski functional of 

Fi +p F2 is proportional to the p-mean of the Minkowski functionals of Fi 
and F2. See p5] for a comprehensive study of p- means of real numbers. 

A convex set K in M"^ is called a convex body if K is compact and has 
non-empty interior. We usually use the letter F for star bodies (which are 
not necessarily convex) and K for convex bodies. 

The support function of a bounded set A in M'^ is defined by 

h{A,u) =snp{{x,u) : X e A} , ueR'^. (2.2) 

Clearly, h{A, u) coincides with the support function of the convex hull of A. 
Any function / : M'^ ^ M, which is sublinear, i.e. f{cx) = cf{x) for all c > 
and f{x + y) < f{x) + f{y) for all x, ?/ G M.'^, is a support function of a convex 
compact set, see |1H1 Th. 1.7.1]. 
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The polar set to a convex body K is defined by 



K* = {u: h{K,u) < 1}. (2.3) 

The same definition apphes if K is not necessarily convex. If K is convex, 
then F = K* is also convex and 

\\u\\F = h{K,u) , ugM'^, 

i.e. the Minkowski functional of F is the support function of K. 

The Firey p-sum of convex sets Ki and K2 that both contain the origin 
can be defined for p > 1 as the convex set L = Ki +p K2 with the support 
function 

h{L, u) = ih{K,, uY + h{K2, uyy/P , (2.4) 

see [in] and [32]. The Firey sum is closely related to the p-star sum (12. ip . 
since 

for convex Ki and K2 that contain the origin. If p = 1, the Firey sum turns 
into the Minkowski sum defined as 

Ki + K2 = {xi + X2 : Xi e Ki, X2 E K2} . 

Then 

h{Ki + K2, u) = h{Ki, u) + h{K2, u) , ueW^. 

Further (without subscript) denotes the Euclidean norm of x G W^. 
The (Euclidean) norm of a set K is defined as \\K\\ = sup{||m|| : u G K}. 
By Volfi{K) or \K\ we denote the d- dimensional Lebesgue measure of K. The 
volume of the unit ^2-ball (i.e. Euclidean ball) B in is denoted by 



i^d = \B\ 



r(i , , 



where F is the Gamma function. The same expression for is used also for 
all real d > 0. 

A random closed set in M'^ is a random element in the space of closed 
sets equipped with the Fell topology and the corresponding Borel a-algebra, 
see [12]. A random closed set X is said to be compact if X has a.s. com- 
pact realisations. If X is a random compact set in R"^ such that ||X|| is 
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integrable, then E h{X, u) is the support function of a convex compact set 
called the (selection or Aumann) expectation of X and denoted by EX, see 
Sec. 2.1]. If X is a simple random convex set, i.e. X takes only a finite 
number of convex compact values Ki, . . . , Kn with probabilities pi, . . . ,Pnj 
then EX = piKi + ■ ■ ■ + PnKn- The expectation of a general X can be 
obtained by approximating X with simple random sets, see Th. 2.1.21]. 

By applying fl2.4p to simple random sets it is possible to define the Firey 
p- expectation EpX, p > 1, of a random compact set X such that G X 
almost surely and E ]|X||^ < oo. In particular, 

h{EpX,u) = (E[/i(X,u)P])'/^ 

is the p-mean of h{X, u) for p > 1, see also ^11] . 

Sets that appear as finite Minkowski sums of segments are called zono- 
topes. Zonoids are limits of zonotopes in the Hausdorff metric, i.e. they 
can be represented as expectations of random segments. By changing the 
Minkowski sum to the Firey p-sum with p>l one obtains Lp-zonoids, which 
appear as limits for Firey p-sums of centred segments. This generalisation 
in the geometric context has been first mentioned in [T7] and has been thor- 
oughly investigated in [36]. It should be noted that Lp-zonoids are exactly 
those sets that appear as polar sets to the unit balls in spaces isometric to a 
(i-dimensional subspace of Lp{[0, 1]), see 

If X is a random closed set with almost surely star-shaped realisations 
(i.e. random star-shaped set), then the p-star expectation of X is the star- 
shaped set F = E* X whose Minkowski functional is given by 

\\u\\F={Y.\\urxf'^. 

This expectation defines a star-shaped set for all p 7^ 0. Note however that 
F is not necessarily a star body, since \\u\\f may be infinite. If X is a random 
convex body that contains the origin and satisfies E ]|X||p < 00, then 

EpX = E;(X*), p>l. (2.5) 
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3 Star bodies associated with SaS distribu- 
tions 



A random vector ,^ G is called symmetric a-stahle (notation SaS) if ^ 
coincides in distribution with — ^ and, for all a, 6 > 0, 

+ 6^/"^ = {a + hf/'^i, 

where ^1,^2 are independent copies of ^, and = denotes equality in distri- 
bution. The value of a is called the characteristic exponent of ^. It is well 
known that ^ is normally distributed if and only if a = 2. 

Theorem 3.1 (see Th. 2.4.3 HZ]). A random vector ^ is SaS with a e (0, 2) 
if and only if there exists a unique symmetric finite measure a on the unit 
sphere S in such that the characteristic function of ^ is given by 

</?5(m) = Ee'<«'"> =exp|-^|(M,z)|V(rf2)| . (3.1) 

The measure a is called the spectral measure of ^. Representation (13.11) 
holds also for a = 2, although the spectral measure is not necessarily unique 
in this case. Although the sphere S can be defined with respect to any chosen 
(reference) norm in M.'^, in this paper we only use the Euclidean reference 
norm, so that S = is the Euclidean sphere in M.'^. 

The expression in the exponential in the right hand side of (13. ip is an 
even homogeneous (of order a) function of u and so defines the Minkowski 
functional of a centred star body F as 

lkllF= [ |(M,z)|V(dz), (3.2) 

SO that 

(^g(u) = e""""^ , ueW^. (3.3) 

The star body F is called the associated star body of ^. Since is finite, 
F always contains a neighbourhood of the origin. If a is not concentrated 
on a great sub-sphere of S"^"^ (i.e. the intersection of S'^"^ with a. {d — 1)- 
dimensional subspace), then ^ is called full- dimensional. In this case ||'u||i;' > 
for all M 7^ 0. 
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The right-hand side of fl3.2l) is called the a-cosine transform of cr, which 
is studied also for all a > —1, see [IB], where it is shown that a is unique if 
a is not an even integer. Although a is not unique for a = 2, the star body 
associated with the normal law is a unique ellipsoid. 

In Section m we see that F is convex if a G [1,2]. Example 14.61 describes 
SaS laws with a < 1 whose associated star bodies are convex. The following 
simple examples deal with general a G (0, 2]. 

Example 3.2 (Complete independence). If ^ is SaS with independent com- 
ponents, then 

i.e. the star body F associated with ^ is ^Q,-ball B^. This star body is not 
convex if a < 1. 

Example 3.3 (Complete dependence). If ^ = (^i,...,.^i) for SaS random 
variable .^i, then 

(p^{u) = exp |- } . 

Note that ^ is not full-dimensional, so that Pf{u) = \\u\\p^ = EmjI is the 
radial function of an unbounded star body F. The corresponding polar set 
F* is the segment with end-points ±(1, . . . , 1). 

The associated star body of SaS random vector ^ can be obtained as 
F = c~^/° E* Yri, i.e. F is the star expectation of 

Y^ = {x: \{x,r])\<l}, 

where c is the total mass of a and rj is distributed according to c^^cr. For this, 
it suffices to note that \{u, z) \ = \\u\\y^ where is the polar set to [—z, z]. It 
is often useful to redefine the spectral measure a to be a probability measure 
on the whole R'^. Then no constant c is needed, so that F = E* F^, where rj 
is distributed in M.'^ according to a. Note that it is always possible to extend 
0" to be a square integrable on R'^, so that the integral (13. 2p exists for all 
a e (0,2]. 

A centred convex body F in R-^ is called an Lp-ball if it is the unit ball of 
a (i-dimensional subspace of Lp([0, 1]). Denote by Cp the family of Lp-balls. 
It is known that F G £p if and only if 

\Hfp= I \{u,z)Y^l{dz) (3.4) 
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for a finite measure fi on S'^"^, see [TSl Lemma 4.8] for p > 1 and [SBJ 
Lemma 6.4] for general p > 0. Note that fl3.4p is called the Blaschke-Levy 
representation of the norm, which is discussed in detail in [25]. It should 
be noted that exp{ — is positive definite for p e (0,2] if and only if 
F G Cp, see |23l [2l] for a survey of related results. By comparing (13.41) with 
Theorem 13.11 and symmetrising, if necessary, the measure fi, we see that 
is exactly the family of associated star bodies of SaS laws. In the following 
we often switch between the letters a and p, since the former is common in 
the literature on stable laws, while the latter is typical in convex geometry 
and functional analysis. 

It is shown in Cor. 6.7] that if F is an Lp-ball with p G (0, 2], then F is 
also an L^-ball for each r G (0,p). It is instructive to provide a probabilistic 
proof of this fact. 

Theorem 3.4. If F is an Lp-ball for p G (0, 2], then F is also an Lr-ball for 
all r G {0,p]- 

Proof. Consider an SaS random vector ^ with a = p and the associated star 
body F. Let C be a non-negative stable random variable with f3 G (0,1). 
Then the characteristic function of ^' = C^'^"^ is given by 

Egi(«',«) = ^-\\u\\f _ 

Thus, F is the associated star body of the symmetric stable ^' with the 
characteristic exponent af3, so that F is an Lr-hall for r = a(3 = p(3 < p. 
Note that ^ and ^' share the same associated star body F. □ 

4 Zonoids and SaS laws with a G [1,2] 

If a G [1,2], it is possible to arrive at a dual interpretation of the charac- 
teristic function (13.11) by noticing that \{u,x)\ is the support function of the 
segment [— x,x], so that 

/ \{u,z)\''a{dz)= [ h{[-z,z],u)''a{dz) = h{K,u)'' , (4.1) 

where 
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is the rescaled Firey a-expectation of the random set X = [—rj, rj] and t] is 
a random vector with values in S'^"^ distributed according to the normahsed 
spectral measure a. The Minkowski inequality implies that h{K,u) is indeed 
a support function of a convex set. If a = 1, then K is called a zonoid, 
see im Sec. 3.5]. Note that representation (14.11) appears already in ^ and 
[2U] in view of its relationship to stable distributions and negative definite 
functions on one hand and Lp-balls on the other one. 

Definition 4.1. Let a be a finite measure on E>'^^^. A convex set K in R'^ 
is called Lp-zonoid with p > 1 and spectral measure a if K = c^^PEp[—r],r]], 
where c is the total mass of a and 1] is distributed according to c~^a. 

It is obvious that Li-zonoids are conventional zonoids. If a is a p- 
integrable probability measure on M'^, then the Lp-zonoid can be defined 
as Ep[— ?7,?7] where rj has distribution a. The following result now becomes 
an easy corollary from Theorem 13. 1[ 

Theorem 4.2. A random vector ^ is SaS with a G [1, 2] if and only if there 
exists a unique centred L^-zonoid K such that the characteristic function of 
^ is given by 

(^^(^z) = e-'^^^'")" , ueR''. (4.2) 

The set K from Theorem 14.21 is said to be the associated zonoid of ^. 
The corresponding polar set F = K* is convex and becomes the associated 
star body of ^. It is well known that all centred convex compact sets on the 
plane are Li-zonoids (i.e. classical zonoids), while this no longer holds in 
dimensions 3 and more. It follows immediately from Theorem 13.41 that the 
family of Lp-zonoids becomes richer if p G [1,2] decreases. 

The following result provides a further interpretation of the well-known 
fact saying that the exponentials of support functions of zonoids are positive 
definite, see Pl p. 194]. 

Corollary 4.3. The function ip{u) = e-'^^^'")", u G M'^, with a G [1,2] and 
a centred convex body K G M.'^ is a characteristic function if and only if K 
is La-zonoid. In this case (p is necessarily the characteristic function of SaS 
random vector. 

A measure a on S'^^^ is called isotropic if the L2-zonoid with spectral 
measure a is a centred Euclidean ball, see [36l 1^ . In other words, if an 
isotropic a is taken as the spectral measure of a Gaussian random vector. 
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then this Gaussian vector has i.i.d. coordinates. The two most common 
examples are the uniform measure on S'^^^ and the cross measure having 
atoms of equal weights at ±ej for the canonical basis ei, . . . , e^. Note that 
the isotropy of o does not mean that the corresponding SaS vectors (with 
a not necessarily equal 2) has a Euclidean ball as its associated star body. 

Example 4.4 (Independent/completely dependent components). The compo- 
nents of SaS vector ^ with a G [1,2] are independent if and only if its 
associated zonoid i^' is a rescaled f^-ball, i.e. 

K = {{aiXi, ttdXd) : x G B'^} 

for ai,...,ad G M. If some of the a^'s vanish, then ^ is no longer full- 
dimensional. Thus, an £g-ball is L,,-zonoid for all r G [l,p] with p being 
reciprocal to q. 

Furthermore, ^ = (ai^i, . . . , a^^i) for a = (ai, . . . , a^) G M*^ and so has 
completely dependent components if and only if K is the segment with end- 
points ±a. In this case ^ is not full-dimensional for each a and d > 2. 

Example 4.5 (Ellipsoids and sub-Gaussian laws). The family of full-dimensional 
L2-zonoids is the family of centred ellipsoids in W^, that also correspond 
uniquely to non-degenerate Gaussian laws on Mf^. Thus ellipsoids are also 
Lp-zonoids for any p G [1,2]. Since polar sets to ellipsoids are again ellip- 
soids, ellipsoids are also Lp-balls for each p G (0, 2]. Ellipsoids do not have a 
unique spectral measure for a = 2. However, if an ellipsoid is represented as 
an Lp-zonoid with p G [1,2) or an Lp-ball with p G (0,2), then its spectral 
measure is unique. The corresponding SaS random vector is said to have a 
sub-Gaussian distribution, see [171 Sec. 2.5]. 

An elliptical norm is determined by a positive definite symmetric matrix 
C, so that \\u\\e = {Cu,u) for the corresponding centred ellipsoid E. A 
simple quadratic optimisation argument yields that 

h{E,u) = yJ{C~H,u), 

see, e.g., pi]. 

Example 4.6 (Sub-stable laws). The distribution of from the proof of The- 
orem [331 is called sub-stable. If ^ is SaS with a G [1,2) and the associated 
zonoid i^T, then ^' is stable with the characteristic exponent a' = af3 and 
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In this case the star body associated with ^' is convex and is equal to the 
polar set to K. In particularly, this holds for all sub-Gaussian distributions 
whose associated star bodies are ellipsoids for each a G (0, 2). 

Theorem 4.7. Each Lp-zonoid with p > 1 and spectral measure which is not 
concentrated on a great sub-sphere of S'^^^ is strictly convex, i.e. its support 
function is differentiable at every point. 

Proof, li p > 1, then | {u, v) \^ is a differentiable function of m, so its integral is 
also differentiable. Since a is full-dimensional, the integral with respect to a 
does not vanish, so that its ^th power is also differentiable. The equivalence of 
strict convexity and differentiability properties is explained in [IBl Cor. 1.7.3]. 

□ 

The strict convexity of K means that for each m G M'^ the support set 

T{K,u) = {yeK: {y,u) = h{K,u)} (4.3) 

is a singleton {x} and the gradient of h{K,u) equals x. Theorem 14. 71 implies 
that polytopes cannot be Lp-zonoids for p > 1, so that the approximation 
by polytopes (often used in the studies of zonoids) is no longer useful for 
Lp-zonoids with p > 1. 

5 Symmetric stable densities 
5.1 Value of the density at the origin 

Consider SaS random vector ^ with a G (0, 2] and the characteristic function 
given by (13. 3p . It is useful to interpret this characteristic function as 

^^{u) = e^ll"!'^ = P{C > \\u\\f} = E Il^>||„,|^ = E l^^^p , (5.1) 

where C is a non-negative random variable with P{C > x} = e~^" for a; > 0, 
so that 

E = r(l + A/a) , A > -a . (5.2) 
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The inversion formula for the Fourier transform yields the following expres- 
sion for the probability density function / of ^ 



= E / e*<"'">dM. (5.3) 

Note that we have used the fact that F is centred. Since / is the expectation 
of the characteristic function of the uniform law on (F, the bounds on this 
characteristic function (see, e.g., [291 Th. 1]) can be used to derive bounds 
for /. 

By substituting x = in (15. 3p we obtain 

Recall that the volumes of F and its polar set K = F* (in case a > 1) are 
related by the Blaschke-Santalo inequality as 

1^1 ■ 1^1 < 4 

with the equality if and only if F is an ellipsoid, i.e. ^ is sub-Gaussian. 

If the spectral measure a is isotropic with the L2-zonoid being the unit 
Euclidean ball and a > 1, then it is possible to apply the results from |36] in 
order to bound the volume of F as 

uJd{'^)/ca < \F\ < ujdia) , (5.5) 

where 

' r(i + • ' r(i + 1) r(fe) ' 

If a e [1, 2), then the equality on the left in (15. 5p is achieved if a is a suitably 
normalised Lebesgue measure on E>'^^^, while the equality on the right holds 
if a is concentrated on ±ei, . . . , ie^. 

5.2 Derivatives at the origin 

Since (p^{u) multiplied by a product of powers of the coordinates of u is 
integrable, representation (15.31) implies that / is infinitely differentiable. Its 
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derivatives at the origin are given by 

Q2m j 



(2vr) 



dxl' ■ ■ ■ dx'/ 



x=0 



a 



where 2m = ki + - ■ ■ + kd. The central symmetry of F imphes that the partial 
derivatives of odd orders vanish. By combining these partial derivatives (with 
m = 1) we arrive at the following expression 



w, 



_d_ 

dxi 



f 



x=0 



(27r 



a 



{w, v)'^dv . 



where w = {wi, . . . ,Wd)- The integral in the right-hand side can be written 
as II w III; where E is an ellipsoid in M*^ called (for a convex F) the Binet 
ellipsoid of F. This elhpsoid is homothetic to the Legendre ellipsoid of F, 
which shares the moments of inertia with F, see j3T|. Results from ^41j can 
be used in order to bound the integral of {w,v)'^ over F. 
Note that 



E 



dxf 



x=0 



(27r 



a 



i=l 



i=l 

where {wlxl + ■ ■ ■ + w'jx'^) defines an elliptic norm of x with the unit ball 
being the centred ellipsoid E with semi-axes iff ^, . . . , W^^. Corollary 2.2a of 
yields that 



[w'^xl + ■ ■ ■ + w'^x'^)dx > 



d 



d + 2' 



'1 [Wi ■■■Wd, 



2/d -2/d 



K 



and so provides an upper bound for the weighted sum of the second deriva- 
tives of the density / of ^ at the origin as 



>f 2d'f 



i=l 



< 



x=0 



And r(i + ^)r(i + f ) 

d + 2 1(1 + ^)1+2/^ 



d\2/d 



with the equality attained if F is a dilate of the ellipsoid E^, i.e. for the 
corresponding sub-Gaussian law. 
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5.3 Expectation of integrable functions 

Integrating (15.31) leads to the following expression 

1 



(2vr)' 



E 



g{-v)dv 



(5.6) 



where g is the Fourier transform of an integrable function g. li g is the 
Fourier transform of a measure /i, then Eg{^) = Fjfi{(F). For example, 
Eexp{ — equals the expected standard Gaussian content of (F. If g 
is the indicator of the Euclidean ball Br of radius r centred at the origin, 
then 

giu) = i2rn/\\u\\r/'j,/,{r\\u\\), 
where Jd/2 is the Bessel function. Therefore 



P{lleil<r}=(-' 



It is also possible to choose g{^) to be the product of functions of indi- 
vidual coordinates of ^, i.e. 



9i0 = ll9i(^^ 



i=l 



For instance, if gi{xi) = li_a,^a,]{xi), i = l,...,d, then 



E 



i=l 



P{^ e xti[-ai,a,]} = 71- 



Ell'^^^^dv. 



i=l 



Vi 



The same argument with the Laplace density gi{xi 

d 



Eexp{- J]A,|^,|} =7r-'^^E 



n 



dv . 



Note that in all these cases the dependency structure is expressed by the set 
F which determines the integration domain, while the value of a influences 
the integrand which is the expectation of a certain function of (. 
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5.4 Renyi entropy and related integrals 

Another instance of (15.61) appears if g is itself the density of Sa'S law with 
associated star body F'. Then 

Eg{0 = ^E\CFnCF'\, (5.7) 
where P{C' > x} = e'^" and (' is independent of (. 

Theorem 5.1. IfC, is SaS with associated star body F and a G (0, 2], then, 
for all c ^ 0, the density f of ^ satisfies 



f{cx)f{x)dx={l + n-"y{Q). 
Proof. Apply (15.71) with g{x) = c'^f{cx) and a = a', so that F' = cF. Then 
EgiO = c' [ f{x)f{cx)dx = E(min(C,cC'))' 



r(i + 4) 



{27rY (1 + c-'^y/" 

Then note that |F| is related to /(O) by flO) . □ 
By choosing c = 1 we see that the density of each SaS law satisfies 

fixfdx = 2-'^/"/(0) . 

The left-hand side can be recognised as the inverse to the 2-Renyi entropy 
power of C,- 



5.5 Probability metric and distance to sub-Gaussian 
law 

Some useful probability metrics are defined using logarithms of characteristic 
functions of random vectors. Extending the definition of the distance between 
two random variables from [521 Ex. 1.1.15] for the multivariate case, it is 
possible to calculate the distance between two SaS vectors ^' and ^" with 
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the same characteristic exponent a G [1,2] and the associated zonoids K\ 
and K2 as 

m^(e',e") = sup{||M||-°|logEe^<"'«'> -logEe*<"'«">| : u ^W^} 
= sup{|/i(/s:i,M)" - /i(K2,m)"| : u G S'^-^}. 

If a = 1, the right-hand side becomes the Hausdorff distance between Ki 
and K2. 

Theorem 5.2. For each SaS vector^ with a G [1, 2) in Mf^ and the associated 
zonoid K there exists a sub-Gaussian SaS vector r] such that ma{^,ri) < 

(c/°/2-l)||if||". 

Proof. For each centred convex body K in M.'^ there exists a centred elhp- 
soid E (called the John ellipsoid) such that E C K C y/dE, see e.g. [T2| 
Th. 4.2.12]. Then it suffices to note that 

\h{K, m)" - h{E, < m)° - h{E, < M)°(ci"/2 _ 1) 

and use the fact that h{E,u) < \\E\\ < \\K\\ for all u G E>'^'\ □ 

Results from Section [6] can be used to related moments of ^ and 77. For 
instance, Theorem l6.ll imphes that E ||^||^/ E ||?7||^ G [l,d^/'^]. 

6 Homogeneous functions of SaS laws 
6.1 Moments of the norm 

If (7 is a homogeneous function, i.e. g{cx) = c^g{x) for all x G and c > 0, 
and so is not integrable over W^, then one can interpret its Fourier transform 
using generalised functions. We refer to [U] for the thorough account of 
generalised functions and their Fourier transforms. The left-hand side of 
(15.61) for not necessarily integrable g can be interpreted as the action of g on 
the density / of ^ (denoted {g, /)), while the right-hand side as the action of 
the Fourier transform ^ of on ip^, i.e. 

is Parseval's identity. Since cp^ given by (15.11) is not necessarily infinitely 
differentiable, the action of a generalised function on it should be interpreted 
as limits if the action of g does not involve differentiation. 
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Theorem 6.1. If ^ is SaS and A G {—d,a), then 

E|Kf-5;r(l±^)g^XjHI^<^- (6.1) 

Proof. Consider (15. 6p for g{x) = \\x\\^ = r^. Using the expression for the 
Fourier transform of g (see \14t, Sec. II. 3. 3]) one arrives at 

E|l«ll' = ^fp|y('-'-''.EiicF), 

where (^r~^~'^, ip) denotes the action of the generalised function r~^~'^ on the 
test function ip. If < A < a, then it is possible to use the regularisation for 
r-^-'^ (see [H Sec. 1.3.9]) to obtain 

(r-^-^EI^^) = EiC')cu, j^^ t-'-\Spit) - l)dt, 

where Ud is the surface area of the unit sphere in M'^ and Spit) is the ratio 
of the surface areas of S{t) fl F and the sphere S{t) of radius t. Then 

rco p poo 

Ud t-^-\l- SFit))dt= / t-^-\l - l[uteF)dtdu 

Jo is'^-i Jo 



oo 

-A-1 



r-'^'dtdu , 



1-1 

If 



which, together with the expression (15. 2p for the moment of (, proves (16.11) 
for A > 0. 

If A G (— d, 0), then no regularisation is needed, so that 

Then (16.11) is obtained by passing to polar coordinates and using the fact 
that r(l - A/2) = (-A/2)r(-A/2). A direct check shows that holds 
also for A = 0. □ 

Remark 6.2. An alternative proof of Theorem 16.11 can be carried over using 
the plane-wave expansion of the Euclidean norm 

I Y{—) f 
=2^^(^f(%X.-J^"'"^"'^"' 
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see [m Sec. 3.10] and the expression for the moments of see Theo- 

rem EUSl 

Example 6.3 (Isotropic law). Assume that ^ is isotropic, i.e. \\u\\p = <7\\u\\ for 
all u and F = B„-i is the ball of radius cr"^. Then (16.11) and the expression 
for the surface area uod of the unit sphere imply that 

E lief = (2a)^ ^-1-^ ^ ff, AG(-ci,a), 

iiMi V ) -p^d^ r(i - I) ^ ' 

which is a well-known formula, see, e.g., [HH Eq. (7.5.9)]. 

If it is difficult to integrate ||m||p' over the unit sphere, it is possible to 
use trivial bounds < \\u\\f < t~^, where R and r are the radii of the 
circumscribed and the inscribed balls to F. The following lower bound is 
sharper for A > 0. 

Corollary 6.4. In the setting of Theoreni \6.1\ with A G (0, a), we have 

with the equahty if and only if F is a Euchdean ball. 
Proof. The expression 

V.x{L,F) = l [ \\u\\-L^-'\\u\\Uu 

is called the dual mixed volume of star bodies K and L (note that the original 
definition [33] is written for the radial functions of L and F). Now it suffices 
to apply the dual mixed volume inequality (see [33] and [351 (2-4)]) 

V_x{L,Fy > |L|'^+^|F|-^ 

with L being the unit Euclidean ball. □ 

Note that the right-hand side of (16. 2p equals E||?7||^, where t] is an 
isotropic SaS random vector with the associated star body being the Eu- 
clidean ball of the same volume as F. 
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Example 6.5 (Multivariate normal and sub-Gaussian distributions). If ^ has 
a multivariate normal distribution with covariance matrix C, then F is the 
ellipsoid E with = \ {Cu,u) and (16.11) implies 

E lief = -^r(^) / {Cu, u)'/'du (6.3) 

for A G (— c?, 2). By passing to the limit, the formula holds also for A = 2. 
The integral retains its value for ^ having a sub-Gaussian distribution with 
the same associated star body F. Thus, the ratio of the moments of the 
norm for a normal vector and the corresponding SaS sub-Gaussian vector 
depends only on a, dimension and the order of the moment. 

If cr^, . . . , 0"^ are the eigenvalues of C, then F = E has semi-axes v^/cj, 
i = 1, . . . , d, whence Kd/\F\ equals 2"^^/^ Yl Corollary 16.41 yields that 

E|,f,2-||)£|i-|(n.)^'^ (6,) 

with the equality if and only if cxi = • • • = 0"^. In particular, if C is diagonal, 
then 

with the equality if and only if E l^f does not depend on i. For this, we 
have used (16.41) and the fact that 

Example 6.6 {SaS vectors with i.i.d. components). Let ^ be SaS with the 
associated star body being ^Q-ball with a G (0,2], so that its coordinates 
^1, . . . , are i.i.d. SaS random variables, see Example 13. 2[ The formula for 
the volume of the £p-ball from [15| p. 11] and (16.21) imply that 



Y(d±\\ p/'l I d\\/d 

E(4i + --- + 4dJ >^ r(^+^) f(^ r(i I '^^ 

for all A G (0, a). For the opposite inequality note that the largest Euchdean 
ball inscribed in Bt has radius d"^. whence 
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for A G (0, a). The both inequahties turn into an equahty for a = 2 and any 
A G (0, 2] and yield the well-known moments of the chi-square distribution 
with d degrees of freedom. 

Using bounds for the average values of norms on the unit sphere from 
it is possible to relate moments of different orders. 



Corollary 6.7. Let ^ be SaS with a G (1, 2] and the associated star body 
F. Let b be the radius of the largest centred Euchdean ball inscribed in F. 
Then for all A G [1, a) 

where Ci and C2 are absolute constants, 

r(— )r(i - -) y§d-i 

and 

2^-ir(^)r(i-A) 



7r'^/2r(l - 1) 

Using [26l Lemma 3.6] for the Fourier transform of the power of the ip- 
norm ||a;||p it is possible to arrive at the following expression 

E mi = ^ ^^^^ I [ s'^'-' n lpisv^)dsdv , 

where is the Fourier transform of the function e~'^'^, x G M. It is valid for 
A G {—d,0) and for A G {0,mm{a,dp)) with non-integer X/p. 

Although II ^11'*' is not necessarily integrable, the integral in the right-hand 
side of (16. ip is well defined for all A > 0. The following result describes the 
limiting behaviour of the A-moment of ||^|| as A t a. 

Corollary 6.8. If ^ is SaS with a G (0, 2) and spectral measure a, then 

Eiier ^ 2^ r(^)r(^)r(f) 
AT.r(i-A) y^r(i-f)r(^)r(^) 
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Proof. It suffices to refer to fl6.1l) together with 



MFdu= / \{u,y)\Mdy)du= / \{u,y)rdu \ aidy) , 

and use the fact that 

/ \(u,y)rdu= , , B(^^—^,-)\\y\r, 

where B is the beta-function. □ 

The intersection body of a centred star body L is the star body IL such 
that 

\\u\\jI =Vo\d-iiLnu^) , ueS'^-^. 

For C e M*^, define = U\\Vo\d-i{L f] ^^), where is the {d - 1)- 

dimensional subspace orthogonal to ^. 

Theorem 6.9. If ^ is SaS with associated star body F and d >2, then 

Proof. It is known [26, p. 72] that the Fourier transform of g{x) = ||a;||^^ for 
a star body L is given by (27r)'^||u||^'^^^/ {n{d — 1)). Thus, 

If F = L, the integral becomes the polar coordinate representation of the 
volume of F. □ 

Similarly to Theorem 16.91 and using [26], Th. 4.6] it is possible to deduce 
that 

where I^L is the fc-intersection body of L, so that this moment is proportional 
to the volume of F is L = F. Note that these intersection bodies are defined 
from Volfc(/fcL fl H^) = Voln-kiL fl H) for each (n — A;)-dimensional subspace 
H, which differs by a factor of 2 from the definition of IL for k = 1. 
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6.2 Mixed moments 

The following result deals with joint moments of the coordinates of ^. For a 
function g{xi, . . . , Xd) and j = 1, . . . ,d denote 

^j9{x) = g{x) - g{x\j) , 
where x\j is x with the jth coordinate replaced by zero. 

Theorem 6.10. If ^ is SaS and Ai, . . . , are positive numbers with A = 
^ Aj < a, then 



i=l " 2 



X / \u,\-^'-'■■■\ud\-^''-\A^■■■AdMu))du. (6.5) 

Proof. The result follows from the formula for the Fourier transform of as 
-2sin(A7r/2)r(A+l)|M|-^-^ (see [H Sec. II.2.3]) and the fact that the Fourier 
transform of the direct product is the direct product of Fourier trans- 

forms, see [m Sec. II. 3. 2]. The expression Ai ■ ■ ■ 1If(m) appears as a result 
of the regularisation procedure, see [HI Sec. 1.3.2]. Finally, one needs the 
expression for the (— A)th moment of ( from (15. 2p and the fact that 

1 . Att A2^^^ r(^) 

— sm— r(A + lj 



TT 2 ^ ' ^ r(i-f) 



□ 



Since Ai ■ ■ ■ A^^ 'Kp{u) vanishes in a neighbourhood of the origin, the in- 
tegral in (16. 5p is well defined. If c? = 2, then (16. 5p turns into 

oA-2 \ ^ X Fl'Ai+i'l 



i=l 

X / \ui\-^'-^\u2\-^'~UlF{Ul,U2)-lF{0,U2)-lF{uuO) + l)dUidU2. 

The signed power of a real number t is defined by 

= |t|^sign(t) , (6.6) 

where sign(t) is the sign of t. 
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Theorem 6.11. If ^ is SaS in for an even d and Ai, . . . , are non- 
negative numbers, with none of them being 1 and such that A = ^ Aj < a, 
then 



X / u['^'''^ ■■■uif^'-'^du, (6.7) 
Jf 



where the integral is understood as its principal value, i.e. the limit of the 
integral over F \ sB as e — >■ . The mixed moments vanish if d is odd. 

Proof. Use the formula 2ir(A + 1) cos(A7r/2)M'>~^~^^ for the Fourier transform 
of the function x^^^ with a non-integer A, see [14^ Sec. II. 2. 3] and identities 
for the Gamma function. □ 

For a centred star body F in M'^ denote 

1(f) 

where the integral is understood as its principal value (note that F contains 
a neighbourhood of the origin in M'^). Note that X{AF) = I{F) for each 
diagonal matrix A with non-vanishing entries and 

J(F)=/ '-^^dv. 

J§k-i Vi-- - Vk 

Corollary 6.12. Jf is SaS random vector in and d is even, then 

7d 



Esign(ei---ed) = (6-8) 



Since the left-hand side of (16.81) does not exceed one in absolute value, we 
obtain an inequality |X(F)| < tc'^ valid for all centred star bodies F C M"^. 
Note that the expectation in (16. 8p does not depend on a. U d = 2, then 

Esign(ei6) = -^2:(F). (6.9) 

Note that in (16. 7p at most one of the Aj's equals 1, since their total sum 
is strictly less than 2. The case of Aj = 1 needs a special treatment, since the 
Fourier transform of x is given by {—2TTi)6'{u), i.e. it acts as {—27ci) times 
the derivative of the test function at zero. 
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Theorem 6.13. Let ^ be SaS with a G (1,2], the associated star body F 
and the associated zonoid K = F*. If d is even and X2, . ■ ■ ,Xd sue non- 
negative numbers such that A = 1 + A2 + ■ ■ ■ + < a, then 

p(A2> p(Ad)^ _ a2^-ii'^ X r(l + ^) 



i=2 ^^2 2! 

I u^^^'''^^ ■■■u^d^''^^\MF~^KT{K,u),ei)du2---dud, (6.10) 

where Ci = (1,0,..., 0), T{K,u) is the support set of K in direction u, 
see ( I4.3I) . and the integral is understood as its principal value. The mixed 
moments vanish if d is odd. 

Proof. Since a > 1, Theorem 14.71 implies that the support function of K 
is differentiable. It is well known (see |48l Th. 1.7.2]) that the directional 
derivative of the support function is given by 

hm h{K,u + vs)-hiK,u) ^ 

slO S 

This formula for v = ei yields that the Fourier transform fi acts on (p{u) = 

(-27rz)e-''(^'"li)"a||u|i||^-^/i(T(ir, u|i), d) , 

where u\i = (0, U2, . . ■ , Ud). The remainder of the proof relies on the formulae 
for Fourier transforms of the signed powers as in Theorem 16. Ill □ 

The following result gives a formula for the probability that SaS vector 
^ takes a value from a polyhedral cone. 

Theorem 6.14. If ^ is SaS with associated star body F, then for each 
invertible matrix A we have 

m=0 {ji,.-,«2m}C{l,...,d} 

where Hi-^^,,,^i^^ is the hyperplane of dimension 2m spanned by the basis 
vectors e^,, . . ■,ei^^. 



28 



Proof. By [HI II. 2. 3 (6)], the Fourier transform of the generalised function 
= 'S.Xj>o is given by iuj^ + '7i6{uj), so that the Fourier transform of I^-gRd 
is the product 

Now it suffices to open the parentheses in the product and use the fact that 
the delta function 6{uk) applied to the indicator of F yields 1. 

Finally, it remains to note that P{^ G AR^} = P{A-^^ G M+} and that 
A~^^ has the associated star body A'^F. □ 

It is easy to see that the result of Theorem 16.141 corresponds to (16. 9p if 
d = 2. In a similar manner it is possible to compute mixed moments of the 
positive parts of the components of ^. 

6.3 Integrals of the density 

The following result expresses the integrals of the density over 1-dimensional 
subspaces of M'^. 

Theorem 6.15. If f is the density of SaS law, then, for each unit vector u, 
I f(^^)dt = 7^:;;^r(l + ^)Af,^ , (6.12) 

where Ap^u = Volrf_i(F fl is the {d — l)-diniensional Lebesgue measure 
of the intersection of F with the subspace orthogonal to u. 

Proof. Using the technique of generalised functions, it is possible to calculate 
the Fourier transform of the function g = 6(^u,x) for a fixed unit vector u as 
{g,ip) = {g,'<P) for any test function ip and its Fourier transform tp, see [Ti] . 
A direct calculation shows that 

{g,ij) = {2ny-' [ ij{tu)dt. 
Jm. 

By applying this expression to the density / and using (15.11) we obtain that 

{g,ip) = EVoh-i{{CF)nu^) = T{1 + ^^)Af,^. 

a 

□ 
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The question, if Ap-^^^ < ^F2,u for convex sets Fi and F2 and all m G S'^^^ 
implies that the volume of Fi is smaller than the volume of F2 is known 
in convex geometry under the name of the Busemann-Petty problem. This 
problem has been recently completely solved (see, e.g. [13\ for the solution 
based on the Fourier analysis) by establishing that the answer is affirmative 
only in dimensions at most 4. The sets F that appear as associated star 
bodies of SaS distributions are Lp-balls and so are intersection bodies, for 
which the Busemann-Petty problem has an affirmative answer in all dimen- 
sions, see [26| Sec. 4.3]. In application to stable distributions this means that 
if two SaS densities /i and /2 with the same characteristic exponent satisfy 



then /i(0) < /2(0). Recall that by (15. 4p the value of the density at the origin 
is proportional to the volume of F. 

It is also possible to consider the intersection of F with a subspace of 
dimension k and obtain that (see also Lemma 3.24]) 



which yields ( 15. 4p for k = d and ( 16.12^ for k = d — 1. For A; = 1 we get 



It is also possible to express the integral of the type f(tu)t'^~^'^~^dt 
by means of the action of the generalised function on the test func- 

tion Ap^uit) = Volrf_i(F n {u^ + tu)). This yields the L^+^-star of ^, see 
[55] . In particular, the Li-star of ^ has the radial function (16.121) and so is 
proportional to the intersection body of F. 

6.4 Scalar products and zonoids of random vectors 

Moments of scalar products of ^ with unit vector u can be easily calcu- 
lated using the Fourier transform of the generalised function \{x,u)\^, see 
[26l Lemma 3.14], or by the explicit calculation of the moments of the SaS 
random variable {^,u). 
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Theorem 6.16. If ^ is SaS and u E E> , then 

Y(^) T(l - -) 

for A G (—1, a). 

The zonoid of an integrable random vector is defined as the expectation 
of the random segment X = [0,^], see [28], Hi]. Representing X as |^ + 
[— 1^, 1^], the expectation of X can be found from 

hiEX,u) = ^{E^,u) + ^E\{^,u)\. 

If ^ is SaS with a G (1,2], then = and Theorem 16. 161 with A = 1 yields 
that 

hiEX,u) = -Til--)\\u\\F, 



TT a 

EX = -r(i - -)K. 

71 a 



so that 



Thus, the zonoid of C, in the sense of Plj coincides with the rescaled asso- 
ciated zonoid of C,- The volume of the zonoid EX is closely related to the 
expectation of a random determinant whose columns are i.i.d. realisations 
of C,- Note also various statistical applications of zonoids of random vectors, 
e.g. for trimming of multivariate observations, see [27]. Furthermore, the 
associated zonoid of ^ can be estimated as the rescaled zonoid of ^, e.g. by 
evaluating the Minkowski average of [0, ^W] for the i.i.d. sample 
In order to recover the spectral measure from the associated zonoid, one has 
to use the inversion formula for the p-cosine transform (see [25]) combined 
with a smoothing operation applied to the support function of K. 



7 Stable laws in 
7.1 Power sums 

Stable laws with all non-negative components (or one-sided laws) are tradi- 
tionally called totally skewed to the right. However, if one considers them 
on the semigroup with addition, these distributions can be also called 
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symmetric stable laws, since this semigroup has the identical involution, see 
[5]. Still we retain the term SaS only for origin-symmetric stable laws in the 
whole space. The Laplace transform of one-sided strictly stable law is given 
by 



E e-<^'«> = exp I - {u, yTaidy) I , ue 



where a G (0, 1) and the spectral measure a on S*^ ^ = §^ ^ fl M^J. is unique. 
It is clearly possible to write 

Eg-(«,f> = e""""^ (7.1) 

for a centred star-shaped (not necessarily convex) set F from (13. 2p with the 
spectral measure obtained by taking all possible reflections of a with respect 
to coordinate planes. 

Below we show how to develop an alternative representation of the Laplace 
exponent using convex sets. For this purpose, it is useful to work with gener- 
alised power sums. For p G (0, cx)), the p-sum of two non-negative numbers 
s and t is defined by 

li p = oo, this operation turns into the maximum of s and t. The p-sum 
X +p y for x,y & M.'^ is defined coordinatewisely as 

X +py = {xi +p yi, . . . , Xd +p yd) ■ 

Random vector ^ in M^^ is strictly stable for p-sums with characteristic 
exponent a 7^ if 

a'/"6+p&'/°6 = + (7.2) 

for all a, 6 > and .^1)^2 being independent copies of ^. The special cases 
correspond to the usual stability for arithmetic sums {p = 1) and max- 
stability {p = 00). The general results from [5] concerning stable distributions 
on abelian semigroups imply that a G (0,p]. It is easy to see that ^ satisfies 
(17.21) with p G (0, 00) if and only if is strictly stable for arithmetic sums 
with the characteristic exponent a' = a /p. Note that a power of a vector is 
always understood coordinatewisely, i.e. ^'^ = (^f, ...,0). The pth signed 
power of a set M C M*^ is defined as 

^{P) = {3,(P> . xeM}, (7.3) 
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where p > and x^^^ is the vector of signed powers of the components of x, 
see fl6.6l) . 

The analytical tools for p-sums rely on the concept of characters on semi- 
groups, see [3]. A character x is a homomorphism between a semigroup and 
the unit complex disk with multiplication operation. The involution oper- 
ation on the semigroup corresponds to the complex conjugation operation 
on characters. The involution is identical if and only if all characters are 
real-valued. In particular, in M+ with (arithmetic) addition (and identical 
involution) the characters are x{^) = e~*^; in M with addition (so that the 
involution is the negation) we set x{^) = ^+ with the coordinatewise 

maximum (and identical involution) the characters are x{^) = ^x<t for t > 0. 

If the characters separate all points, i.e. if Xi ^ x^ implies x{x\) ^ 
for some x and the characters generate the Borel cr-algebra, then the Laplace 
transform x ^ Ex(0 characterises uniquely the distribution of a random 
element ^, see O Th. 5.3]. In special cases one obtains the characteristic 
function, the Laplace transform or the cumulative distribution function. 

In Mi with the p-sum operation the characters are given by 



d 



Xu[,x) = exp 



J2{^^Uiy\ , xeRl, (7.4) 



1=1 



for u G M'^ if p is finite, li p = oo, the characters are 

. . I 1 if XiUi < 1 for alH = 1, . . . , d, . 
Xu{x) = < . (7.5) 

U otherwise , 



for u e Rl, so that Ex«(0 = < with m"^ 



1^1 >--->Wd 



7.2 Li(j9)-zonoids 

Let 

yM = {{yixi, . . . , ydXd) : x E M} (7.6) 
denote the set M C M'^ rescaled by a vector ?/ G M'^ and a set M C M°'. 

Definition 7.1. Let a be a finite measure on S*^^^ with total mass c. Define 
to be a random vector distributed according to c^^cr. The set K = cYiX 
for 

X = 7]B^„ = {{t],v,, . . ..riava) : \\v\\, <l, v e (7.7) 
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with p ^ + q ^ = lforp>lis said to be Li{p)-zonoid with spectral measure 
a. 

Definition 17.11 can be reformulated for a probability measure a on M*^ and 
the corresponding random vector 77. In this case K = E{r]B^) is also called 
the Li(p)-zonoid generated by rj. 

Note that X from (17.71) is a rescaled ^^-ball B^. If p = 00, then X becomes 
a rescaled crosspolytope. More generally, taking the Firey a-expectation 
Eq, X yields an LQ,(p)-zonoid. The conventional and Lp-zonoids are not mem- 
bers of these new families. It is however possible to define a family of sets 
that includes all zonoids introduced so far. 

Definition 7.2. Let M be a centred star-shaped set in M'^ and let 77 be a 
random vector in M'^ with E||r7||^<oofor]9> 1. The Firey p-expectation of 
T]M is called Lp{M)-zonoid. 

If M is the segment with end-points ±(1, . . . , 1), then Definit ion 1 7 . 2 1 yields 
the family of Lp-zonoids. The case of M being simplices of varying dimen- 
sion was considered in [16] . If M is an £g-ball, we arrive at Definition 17.11 
Although Definition 17.21 with a general M may be of geometric interest, we 
do not pursue its study in this paper. 

It is obvious that Li(p)-zonoids are plane-symmetric, i.e. they are sym- 
metric with respect to all coordinate planes. The following proposition 
shows that the Li(p)-zonoid is actually determined by the vector Ir^l = 
• • • , \Vd\) of the absolute values of r] = (r/i, . . . ,rid). This means that 
it suffices to consider only spectral measures on M^. 

Proposition 7.3. If\ri'\ and |?7"| share the same distribution, then the Li{p)- 
zonoids generated by Vj and rj" coincide. 

Proof. It suffices to notice that rj'B^ and T]"Bg coincide in distribution. □ 

Theorem 7.4. A random vector ^ G MJ^ is strictly stable for p-sums with 
a = 1, p G (1, 00] and spectral measure a if and only if 

Ex.(0=e-'^^''•"^ (7.8) 

for an Li{p)-zonoid K with spectral measure V{1 — ^)a, where Xu is given 
by ( I7.4|) ifp is finite and by ( I7.5|) if p = 00. 
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Proof. Assume that p G (l,C)o) and consider a G The general results 

from [51 Sec. 5.3] imply that the Laplace transform of a strictly stable for 
p-sums random vector in Ml with characteristic exponent a is given by 



where 

^{u)= f {l-Xu{x))k{dx), mgM^, (7.9) 

and A is the Levy measure ^. The Levy measure admits the polar decom- 
position as at~'^^'^dta{dy) for x = ty, so that a change of variables in the 
integral yields that 

\ a/p 

Tiu^yA a{dy)- H {1 - e-^)s-f-'ds . 



Thus, 

Y^i^^yA <y{dy). (7.10) 
1=1 ) 

Since the £p-norm of u can be written as (Xl'^D^'''' — hiBg,u), 



where yB^ is defined as in (17. 6p . If a = 1, then 

^{u)=r{l--) [ h{yB'^,u)a{dy), 

i.e. V(m) = h{cEX,u) for u G R'^, where c = cr(§"[r^)r(l - i) and X given 
by (17. 7p with j9~^ + = 1 and ?7 distributed according to the normalised 
a. The case p = oo is considered similarly, see [13] for the special study of 
max- stable laws. □ 

Remark 7.5 (Arithmetic sums). Although the conventional case of arithmetic 
sums (p = 1) is not covered by Theorem 17.41 (17. 8p also holds. Then q = 
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oo, so that dlZD reads X = xf^J-r/^, r^^] for rj G M^. Thus K = EX 
xf=i[-Er/i,Er7i], so that 



i.e. ^ is deterministic. Indeed, one-sided strictly stable laws with a = 1 are 
necessarily degenerated. 

Example 7.6 (Max-zonoids) . Let ^ be a max-stable random vector in M^^, 
whose marginals are unit Frechet, i.e. ^ is stable with respect to the coordi- 
natewise maximum operation and exponent a = 1. The Laplace transform 
of ^ is the cumulative distribution function at the point and 

for an Li(cxo)-zonoid K, i.e. K is the expectation of the randomly rescaled 
crosspolytope in W^. These zonoids (more exactly their intersections with 
M^) have been explored in [43j in view of the studies of max-stable distribu- 
tions, and so are called there max-zonoids. 

Example 7.7 {p = 2). If ^ is strictly stable for p-sums with p = 2 and a = 1, 
then the Laplace transform of ^ is given by 



h{K,u) 



The Li(2)-zonoid K is the selection expectation of a randomly rescaled Eu- 
clidean ball, i.e. the centred plane-symmetric random ellipsoid. 

Example 7.8. In dimension c? = 2 it is possible to calculate the support 
function of X from (17. 7p for u = {ui, U2) € as 



/i(X, u) = sup{Mir/i cos^/" e + U2ri2 sin^''^ 9: 0<e<-}. 

Substituting the value of the critical point tan^ = {u2r]2/uir]iy/^ and notic- 
ing that 'r]i,ri2 > we arrive at 

h{K,u) =Eh{X,uj =E 



[mi?7i + U2r]2 



\p-i 



Example 7.9 (Completely dependent and independent cases). If t] is deter- 
ministic, then the corresponding Li(p)-zonoid is a rescaled £q-ball and the 
coordinates of ^ are completely dependent. 
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A centred parallelepiped x^^^[— aj,aj] is an Li(p)-zonoid for each p > 1. 
To check this, it suffices to take the spectral measure concentrated at the unit 
basis vectors ei, . . . , with masses oi, . . . , a^, so that ^ has independent 
components. For instance, X from (17.71) becomes the segment [— ej,ej] if 
rj = Ci. Therefore, polytopes may be Li(p)-zonoids for p > 1, cf Theorem 14. 71 

Thus, r] equal to one of the basic vectors results in X from fl7.7p being a 
segment, while any 7] from the interior of results in X being a rescaled 
£g-ball. By combining such values of rj it is easy to construct further ex- 
amples of Li(p)-zonoids. For instance, if t] takes the values (2, . . . , 2) and 
(2, 0, ... , 0) with equal probabilities 1/2, then the corresponding Li(p)-zonoid 
is the Minkowski sum of the unit £q-ball and the segment with end-points 
±(1,0,. ..,0). 

7.3 One-sided strictly stable laws 

The construction based on p-sums makes it possible to provide a geometric 
interpretation of strictly stable laws for arithmetic sums on and a G (0, 1]. 

Theorem 7.10. A random vector ^ G is strictly stable (for arithmetic 
sums) with a G (0, 1] if and only if the Laplace transform of ^ is given by 

Ee-«.«> = e-'^^^'"") , ueRi, (7.11) 

where = {u", . . . , u^) and K is Li{a~^)-zonoid called the associated zonoid 
of ^. The spectral measure of K is r{l — a)a, where a is the spectral measure 

Proof. The random vector is strictly stable for p-sums with p = ^, so that 
(17. lip follows from Theorem 17. 4[ If a = 1, the law of ^ is degenerated, see 
Remark 17. 5[ □ 

In particular, if a = |, then K from (17.111) is the expectation of a random 
ellipsoid as in Example 17. 7[ By comparing (17. lip with (17. ip we see that 

\\u"\\f* = so that K* = 

Example 7.11 (One-sided sub-stable laws). Let C, be one-sided stable law 
in with a G (0, 1) and let C be a non-negative stable random variable 
with characteristic exponent P G (0,1). Then ^' = C^^"^ has a sub-stable 
distribution, see Example 14. 6[ It is easily seen that 
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where K and L are the associated zonoids of ^ and ^' respectively, i.e. K is an 
Li(a~^)-zonoid and L is an Li((a/?)~^)-zonoid. Note that \\u\\^p = \\u^\\p{i3) 
for any star body F, where F^^'^ is defined by (17. 3p . Hence h{K,u°')^ = 
h{L,u''^) where L* = {K*Y^\ 

Theorem 7.12. If K is an Li{p)-zonoid for p > 1, then K is Li{r)-zonoid 
for all r > p. 

Proof. We refer to the construction from Example I7.11[ Assume that K is 
a parallelepiped, i.e. the components of ^ are independent. Then L* is the 
/3-power of the crosspolytope K*. Since the crosspolytope K* is the (possibly 
rescaled) £i-ball, its /?-power (K*Y^'^ is the (possibly rescaled) £i/^-ball. Its 
polar L = ((fC*)^^^)* is a (possibly rescaled) £i/(i_/3)-ball. By the construction 
of Example 17. m the £i/(i_^)-ball L is an Li((a/5)~^)-zonoid for all a G (0, 1). 
By setting q = 1/(1 — it is easy to see that £g-ball is Li(r)-zonoid for all 
r > p, where p~^ + = 1. 

Thus, ig-hal\ can be represented as the expectation of rescaled ^^'-balls 
for each r' < q. Since each Li(p)-zonoid is the expectation of the rescaled 
iq-hal\, it can also be expressed as the expectation of rescaled ir'-haW, where 
r' is associated with r, so that it is also an Li(r)-zonoid. □ 

Thus, the family of Li (j9)-zonoids becomes richer if p increases. The 
richest one is the family of Li(oo)-zonoids (or max- zonoids) . In the planar 
case it includes all plane-symmetric convex sets [33], while in the spaces of 
higher dimensions this is no longer the case. 

The following generalisation of Theorem 17.41 treats general strictly stable 
laws for p-sums. 

Theorem 7.13. A random vector ^ G is strictly stable for p-sums with 
p G (0, oo] and the characteristic exponent a <p if and only if 

Ex„(0 = exp{-/i(ir,n")} 

for an Li{p/a)-zonoid K. 

Proof. If ^ is strictly stable for p-sums with a < p and finite p, then is 
Sa'S for arithmetic sums with a' = a/p G (0, 1], so that Theorem 17. lUI yields 
that 

Eexp|-5^(e.M,)^'} =exp{-/i(K,M"'^)} =exp{-/i(K,M")}, 
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where K is an Li(p/a)-zonoid. The case of p = oo follows from the fact that 
ip{u) from (17.91) is written as 

V^(m) = / {uvfxx{uiyi)Ya{dy) = / max (wjyf )a'((i?/) 
for another measure a' . □ 



7.4 Moments of one-sided stable laws 

Similar to the Fourier analysis technique in Section Ej it is possible to use the 
expression for the Laplace transform of one-sided stable law in order to obtain 
more information about its probability density function / and moments. By 
integrating the both parts of 

e-<"'">/(x)rfa; = e-'^^^'"") , m 6 M!^ , (7.12) 

with respect to u with a certain weight, we arrive at various expressions for 
the moments of ^. To start with, integrate the both sides of (I7.12p over the 
ray {tu : t > 0} with weight for a fixed unit vector u and A > — 1. Since 

oo 

^-{x,u)t^x^^ = {x,u)-^-^T{l + A) , 





the integration of the right-hand side yields that 

a 1 (1 -|- Aj 

For instance, if K is the parallelepiped y<f=i[~ai,ai] (which corresponds to 
the independent coordinates of ^) , then for A = we have 



B«,„)-..ra.i,(E.<)-""^(Efe)°) 



-l/a 



Note that we have used the fact that E,^~^ = r(l + l/a)a- ^^'^ . 

If one performs a similar integration with A G (— 1 — a,— 1), it is possible 
to regularise the integral of e"^^'""*^ by subtracting the value of the function 
for t = as 

/ (l-e-<^'">)/(x)dx = l-e-"(^'^"). 
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Then, for /3 G (0, a) 



Clearly, the above expressions for the moments can be obtained by calcu- 
lating the moments of one-sided stable random variable {^,u). Furthermore, 
similar results can be obtained for random vectors which are strictly stable 
for p-sums. The case of p = oo is considered in |43j . 

The multivariate Laplace ordering is introduced in [49] by pointwise or- 
dering of the Laplace transform. Thus, two one-sided strictly stable random 
vectors with the same characteristic exponent are Laplace ordered if and 
only if the corresponding associated zonoids are ordered by inclusion. Ap- 
plications of this ordering for actuarial quantities have been considered in 

8 Geometric interpretations of the spectral 
measure 

8.1 j9-surface area measures and spectral bodies 

Assume that the SaS distribution is full- dimensional, i.e. the spectral mea- 
sure a is not concentrated on a great sub-sphere of S"'"^. The Minkowski 
existence problem [HI Sec. 7.1] establishes that for each finite positive full- 
dimensional even Borel measure on the unit sphere there exists a unique 
centred convex body Q such that S{Q,-) = cr(-). Here S{Q,-) is the sur- 
face area measure of Q which is the unique measure on the unit sphere that 
satisfies 

dV^iQ, L) = lim 1^ ^ = I h{L,u)S{Q,du), 

elO e Jsd-1 

where Vi{Q,L) is called the mixed volume of the convex sets Q and L. We 
refer to [IE] for a detailed presentation of the relevant concepts from convex 
geometry. 

The Lp generalisation of the above concepts has been studied in [321 E3] ■ 
The p-mixed volume Vp{Q,L) of two convex bodies containing the origin is 
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defined by 

-VJQ, L) = lim \Q+p^^'''L\-\Q\ ^ 1 /■ ^(^^ ufSJQ, du) , (8.1) 

P elO e p Jgd-1 

where the p-surface area measure Sp{Q,-) satisfies 

Sp{Q, du) = h{Q, uy-PS{Q, du) . 

The p-Minkowski problem is solved in [321 Th. 3.3] by showing that if a is 
an even positive Borel measure which is not concentrated on a great sub- 
sphere of S"^"^ and p > 1, p d, then there exists a unique centred convex 
body Q, such that Sp{Q, ■) = cr(-). By combining this representation with the 
classical Minkowski problem for p = 1, any full-dimensional spectral measure 
a corresponding to SaS law ^ with a G [1, 2) (a G [1, 2] in dimension d >3) 
can be interpreted as the a-surface area measure of a centred convex body 
Q, i.e. (t(-) = Sa{Q,-)- We call Q the spectral body of ^. By (13.21) . the 
Minkowski functional of the associated star body F (or the support function 
of the associated zonoid K) can be expressed as 

||M||^ = /i(ir,u)"= / \{u,v)\''S^{Q,dv). 

The pth projection body HpQ of Q is defined in [3l] by 

h{UpQ, u)P = — ^ / I {x, u) \PSp{Q, dx) , 

dHdCd-2,p J§d-i 

where 

The set HpQ (or its dilated version) is sometimes denoted by T*_pQ and is 
called the polar centroid body, see, e.g. [37], (4.2)]. The normahsing constant 
guarantees that UpB = B for the unit Euclidean ball. Thus, the associated 
zonoid K satisfies K = {dKdCd-2,aY^°'^pQ- 

The Lp-analogue of the Petty projection inequality proved in [3^ Th. 2] 
estabhshes that 

\Q\id~P)/P.\YV;Q\<K''Jr 

Using the fact that the polar set H* Q is a dilate of the associated star body 
F of ^ and setting p = a, we arrive at the following inequality 

|Qr'+'/" ■ |i^| < (rfQ_2,„)-"/" (8.2) 
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valid for a G [1, 2] with the equahty attained in the sub-Gaussian case. Recall 
that \F\ determines the value of the density of the stable law at the origin 
and provides a bound for the moments of ||^||. In the Gaussian case a = 2 
for (i > 3 and (18.21) reads 



-d/2 



8.2 Spectral star body 

The following interpretation of the spectral measure is useful for SaS laws 
with arbitrary a G (0,2]. Assume that the spectral measure o" of ^ has a 
positive continuous density on S''"^ with respect to the (d — l)-dimensional 
surface area measure, so that 

a[du) = -^—pM'^^'^du 
d + a 

for a star body L, which we call the spectral star body of ^. By passing to 
polar coordinates it is easily seen that 

\Mf= [ \{u,y)ra{dy)= [ \{u,y)rdy. 

The integral in the right-hand side is related to the p-centroid body TpL and 
its polar r*L defined (up to a possibly different normalisation) by 

h(TpL,uY = \\u\\^,^ = — \— I \{u,y)\Pdy, 

see [12], |3l] and [35l (6.1)]. Note that TpL is convex for all p > 1. Thus, the 
associated star body of SaS vector is related to its spectral star body by 

F=(Q,a|L|)-i/" TIL. 

It is proved in [31] that |rpL| > \L\ if p > 1, which implies the Blaschke- 
Santalo inequality \L\ ■ \T*L\ < with equality if and only if L is a centred 
ellipsoid, see also [38]. If p = a > 1, then 
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The same spectral star body L can be used to construct SaS random 
vectors ^oi,l with varying characteristic exponent a. If L is convex, then [151 
Prop. 2.1.1] yields that there exists a universal constant c > such that for 
all u eW^ and p > 1 

\{u,y)\^dy\ <cp [ \{u,y)\dy, 



where Li = \L\^^^'^L has volume 1. Thus 



where Fa^i is the associated star body of C,a,L- By Theorem 16. H this yields 
an inequality between the A-moments of the norm with A G (0, 1) for an SaS 
random vector with a > 1 and the Cauchy random vector with the same 
spectral star body. For the Cauchy distribution we have a = 1, where a 
number of further inequalities for the volumes of projection and polar bodies 
are available, see [T2| H8]. 

If F is the associated star body of an SaS law with convex spectral body 
L, then [131 Lemma 3.1.1] implies that 

For instance, this inequality may be used in order to obtain lower bounds for 
the moments of ||^||. It also means that the associated zonoid of ^ contains 
a dilate of L. 

It is likely that other geometric properties, e.g. curvature, surface area, 
other intrinsic volumes, of the spectral sets and associated star bodies have 
a bearing in view of the studies of SaS laws. 



8.3 Spectral sets for one-sided laws 

Let ^ be a one-sided strictly stable random vector. Although its spectral 
measure a is supported by S+~^, consider its extension on the whole S^"^ in 
a plane-symmetric way. Note that this extension is full-dimensional, so that 
the Minkowski existence problem guarantees that a{-) = S{Q, ■) for a convex 
plane-symmetric body Q, also called the spectral body of ^. 
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Let M be a centred convex set in W^. If 77 is a random vector distributed 
according to the normalised S{Q, ■) having the total mass c, then the expec- 
tation of crjM is given from 

h{cEriM,u)= [ h{yM,u)S{Q,dy) 

h{uM, y)S{Q, dy) = dVi{Q, uM) . 



If ^ is one-sided strictly stable random vector with characteristic exponent a 
and spectral measure a, then choose M to be the £i/(i_Q,)-ball, so that 

Ee-<«'"> =exp{-r(l-a)dyi(Q,u5i/(i_,))}. 

The first Minkowski inequality Vi{K, L) > \K\^'^-^^/'^\L\^/'^ (see HH Th. 6.2.1]) 
together with the formula for the volume of the unit £p-ball (see ^5', p. 11]) 
imply that 

Ee-<«'"> < cxdI ^{'^-(^f \Q\{d-i)/d \ 

^""""^X r(rf(i-«))v<^'^' • 

9 Covariation and regression 
9.1 Bivariate case 

The covariation replaces the concept of covariance for SaS vectors. If ^ = 
('^1,^2) is SaS in with a > 1 and the spectral measure a, then the 
covariation of ^1 on ^2 is defined by 



[6, 6] a = / sisf ^^a{ds) 



see [m Sec. 2.7]. It is mentioned in [HI Sec. 2.7] that the covariation can be 
equivalently defined as 



1 da-{t,M) 



(9.1) 

tl=0,t2 = l 



a dti 

where a(ti, ^2) is the scale parameter of y = ti^i + ^2^2, i-e. 

c^°(ti,t2)=/ \tiSi+t2S2\''a{ds). (9.2 
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Theorem 9.1. If ^ = (^1,^2) is SaS with a G (1,2] and the associated 
zonoid K, then 

[^l,^2]a = XiXr\ (9.3) 

where T{K, (0, 1)) = {(xi,X2)} is the support point of K in direction (0, 1), 
see M) . 

Proof. By Theorem 14. 7[ Lp-zonoids with p > 1 are strictly convex, so that the 
support set T{K, u) is indeed a singleton for each direction u. The right-hand 
side of fl9.2p can be identified as h{K,u)" for K being the associated zonoid 
of ^. The partial derivative in the right-hand side of (19 .ip then becomes the 
directional derivative of h{K,u) in direction (1,0). By [48, Th. 1.7.2] this 
derivative can be expressed as h(T{K, (0, 1)), (1, 0)). Hence 



[^1, 6]. = HK, (0, l)r-'h{T{K, (0, 1)), (1, 0)) = x,x 
It is shown in [47l Lemma 2.7.16] that, for all p G (1, a), 



□ 



(9.4) 



2\a 



Using Theorems 16.101 and 16.131 it is possible to calculate the moments in the 
left-hand side explicitly as 

0-1), _ «2^-V,o P^ r(l + ^) 2x,xl~' 



Eiuri = ^r(2 



rfi±£) rfi - ^) 

Ei6r = 2M^-i-^, 

v^r r(i-f) 

where xi and X2 are the coordinates of T{K, (0, 1)). By dividing these ex- 
pressions we arrive at Xi/x2, which is exactly the right-hand side of (19. 4p . 



9.2 Multivariate case 

The following result provides covariations for random variables that belong 
to a linear span of an SaS random vector. 
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Theorem 9.2. Let ^ be SaS in with a E (1,2]. If u',u" are non-zero 
vectors in M.'^, then 

[{^,u'), = hiK,u'r^'hiTiK,u"),u') . (9.5) 

Proof. The scale parameter of ti{C^,u') + t2{^,u") is h{K,tiu' + t2u"). By 
differentiating its power with respect to ti as in (19.11) . we arrive at (19. 5p . □ 

Theorem 19.21 provides an ahernative reformulation of [17} Lemma 2.7.5]. 
The right-hand side of (19.51) considered a function of u' is the support function 
of the singleton h{K,u")°'''^T{K,u"), and so is additive with respect to u'. 
In particular, if u' = (1, 1, 0) and u" = (0, 0, 1) in R^, it yields the additivity 
of the covariation of SaS random variables with respect to its first argument. 
Similarly, one deduces the additivity of the covariation with respect to the 
sum of independent second arguments. Furthermore, the covariations in the 
left-hand side of (19.51) for all u', u" determine uniquely the associated zonoid 
K. 

Example 9.3 (£p-balls). Assume that K is the unit ^Q-ball, which corresponds 
to SaS vector ^ with i.i.d. components. The support point T{K, u) equals 
the gradient of = h{K,u), see [HI Cor. 1.7.3]. Therefore T{K,u) = 
{\\u\\l-"u^^-'^}. By (inSD, 

9.3 Regression coefficients and linearity conditions 

The covariation is used to build regression models for SaS distributions. By 
dZl Th. 4.1.2], 

E(ei|6) = |Ht^^2 a.s. 
Since [,^2, C2]a = 2^2 T{K, (0, 1)) = {(a;i, X2)}, we obtain 

E(eiie2) = -6 a.s. 

Thus, the regression line is the line passing through the origin and the support 
point T{K, (0, 1)). Therefore, the regression lines are identical for any two 
SaS laws that share the same associated zonoid, e.g. for the Gaussian law 
and its sub-Gaussian variant. 
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It is known [171 Sec. 4.1] that multiple regression is not always linear for 
a e (1,2). The necessary and sufficient conditions for the linearity given in 
[lO] can be reformulated geometrically as follows. 

Consider a convex set K and the one-dimensional subspace spanned by 
X G M°'. The shadow boundary of K in direction x is the set d{K + Hx) (IdK, 
where d denotes the topological boundary in W^, see \5U\ Def. 3.4.7]. 

Theorem 9.4. Let (^i, . . . , C,d) be an SaS random vector with a G (1, 2] and 
the associated zonoid K. Then E(^i|^2, • • • , ^d) hnear in ^2, ■ ■ ■ if ^.nd 
only if the shadow boundary of K in direction ei = (1,0,..., 0) is a subset 
of a {d — 1)- dimensional hyperplane, which does not contain ei. 

Proof. By Theorem 3.1 from [ID], the conditional expectation is linear if and 
only if, for all M2, . . . , Ud, 



By differentiating (14.21) and using [IHl Th. 1.7.2] for the directional derivative 
of the support function, it is easily see that this holds if and only if 



where m|i is m with the first coordinate replaced by zero. Since h{T{K, u\i), Ci) 
is the ith coordinate of the singleton T{K,u\i), (19.61) means that T{K,u\i) 
is orthogonal to a = (1, —02, . . . , — a^) for all u\i = (0, U2, ■ ■ ■ , uj). In other 
words, the shadow boundary of K in direction ei lies in the hyperplane 
orthogonal to a. By the condition, the ffist coordinate of a is not zero, so 
that this hyperplane does not contain ei. □ 

Example 9.5 (Sub-Gaussian laws). If ^ is sub-Gaussian with the norm = 
{Cu, u), then K = E* is an ellipsoid. It is easy to see (see [21]) that T{K, u) 
is the point C~^u/ a/ {C~^u, u). Then the condition of Theorem 19.41 holds 
with a = Cei. 

Corollary 9.6. Let (^1, . . . ,C,d) be an SaS random vector with a G (1,2], 
the associated zonoid K and the spectral measure a. Then E(^i|^2, • • • , ^d) Is 
linear in C,2, ■ ■ ■ ,^d if s^nd only if there exists a G M*^ with non-vanishing first 




d 




(9.6) 



i=2 
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coordinate such that one of the following equivalent conditions holds for all 
u orthogonal to Ci. 



(grad h{K,u),a) = , 



(9.7) 



/ {y,a){y,uy''-'^a{dy) = 0, 



(9.8) 



[(e,a),(e,n)]„ = 0. 



(9.9) 



Proof. Since the support function of K is differentiable by Theorem 14.71 the 
support point of K in direction u is given by the gradient of h{K, u), see [iSf 
Cor. 1.7.3]. This yields, (19. 7p . By differentiating (14.11) . it is easy to see that 



so that (19. 7p is indeed equivalent to (19.81) . Finally, (19.51) imphes that 



The vector ^ has the multiple regression property if, for each linear trans- 
formation A, the multiple regression of the first coordinate of rj = AC, onto the 
remaining coordinates is linear. Note that rj has the associated zonoid AK. 
By Theorem 19. 4[ this happens if and only if the shadow boundary of K in 
each direction is contained in a (c/— l)-dimensional hyperplane. W. Blaschke 
proved in 1916 that for dimension d > 3 this is the case if and only if K is an 
ellipsoid, see [5Ul Th. 3.4.8]. By Theorem 19. 4[ this geometric result translates 
into the multiple regression criterion from pf7'. Prop. 4.1.7]. 

10 Operations with associated sets 

If ^' and C," are independent SaS with associated star bodies Fi and F2, then 
^ = ^' + ^" has the characteristic function 



Thus, ^ has the associated star body F being the radial sum (also called 
a-star sum) of Fi and F2, i.e. F = Fi F2. 





□ 



Ee'<"'«> =exp{-(||M||^^ + \\u\\%J}. 
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Theorem 10.1. If ^' and ^" are independent SaS in M'^ with a G [1, 2] and 
probability densities f^i and f^n respectively, then the probability density of 
^ = ^' + ^" satisfies 

/c(o)-"/'^>/e(o)-"/'^ + /5"(o)-"/'^ 

with the equality if and only if the associated star bodies of ^' and ^" are 
dilates. 

Proof. The resuh follows from (15.41) and the dual Brunn-Minkowski inequal- 
ity for radial sums of star bodies, see |33, Prop. 1.12]. □ 

The following result concerns approximation by sub-Gaussian laws. 

Theorem 10.2. A law is SaS with a G [1, 2] if and only if it can be obtained 
as a weak limit for the sums of independent sub-Gaussian laws with the same 
characteristic exponent. 

Proof. By fTEl Cor. 6.14], each centred convex body F is an Lp-ball with 
p > 1 if and only if can be uniformly approximated for u from the unit 
sphere by finite sums of the form \\u\\^^_^ -|- ■ ■ ■ -|- HmH^;^, where Ei, . . . , 
are centred ellipsoids. The proof is completed by setting p = a and using 
the fact that exp{ — is the characteristic function of a sub-Gaussian 
law. □ 

Consider maximisation of E | (^, m) over m G M'^ for fixed A G (0, a) under 
the constraints {u, fi) = r for some /i G M.'^, r > 0, and {u, (1, . . . , 1)) = 1. 
By Theorem 16. 161 this is equivalent to maximising \\u\\f for u satisfying the 
constraints, i.e. its solution is the direction of the smallest radius-vector 
function for the set F (1 H, where 

H = {ueR'': (n,/i)=r, (m, (1, . . . , 1)) = 1} . 

This corresponds to the idea of portfolio selection studied in [2] for a > 1. 

It is also possible to consider further optimisation problems for the mo- 
ments of the norm of 1] = A^, where A is an invertible linear transform and 
a G [1,2]. The direct computation shows that rj has the associated zonoid 
AK and the associated star body By Theorem 16.11 minimising 

E ll^^ll'*' for A G {—d,a) over A G SL^, corresponds to the minimisation of 
the integral of ||y4^u||^ over u G 
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Consider a special case of this problem for A = 1 and a G (1, 2]. In terms 
of the associated zonoid K = F* , we can equivalently minimise the mean 
width 

w{AK) = 2 / h{AK,u)du, 

over A G SL„. It is shown in [16] that AK has the minimal width position if 
the measure on the unit sphere with density h{AK, ■) is isotropic. 

Taking a subvector of ^ corresponds to a section of the associated star 
body F by the corresponding coordinate subspace. By applying orthogonal 
transformations, we see that the projection of ^ on any subspace H has the 
associated star body {F n H) Q) H^, i.e. the direct sum of F D H and the 
space orthogonal to H. Therefore, the values at the origin of the probability 
density function of the projected ^ are closely related to the intersection body 
of F. 

A bound on the volume of a convex set using volumes of its {d — 1)- 
dimensional sections (see [32] and [121 P- 341]) yields the following inequality 
for the values of the density function at the origin 

where /-j(0) is the density at the origin of the subvector of ^ with the ith 
coordinate excluded. This inequality holds for all SaS laws with convex 
associated star bodies. In the bivariate case, 

r(i + ^) 

where /i and /2 are the marginal densities. Note that the coordinates of 
^ = (^1, can be independent with a convex F only if a G [1, 2]. 

Theorem 10.3. Consider two SaS random vectors ^' and C," of dimensions 
di and d2 and the random vector rj = {^',^") obtained by concatenating ^' 
and Decompose M.'^ into the direct sum of two hnear subspaces Hi and 
H2 of dimensions di and d2 ■ 

Then ^' and ^" are independent if and only if the associated star body F of 
r] (or associated zonoid KifaE [1, 2]) is the a-star sum of Fi = (FnHi) x H2 
and F2 = Hi X (F n H2) (respectively K is the Firey a-sum of the projection 
of K onto Hi and onto H2 if cx G [1, 2]j. 
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Proof. Note that ^' and C 



are independent if and only if 



ll(«l,«2)||^ 



Ml G /fi, U2 e H2. 



It remains to observe that the associated star bodies of subvectors appear 
as the intersections of F with Hi and H2 and the associated zonoids are 



The duahty operation transforms convex bodies into their polar sets. This 
operation does not generally preserve the property of a set being a zonoid 
or Lp-ball. However, it makes sense if applied to the spectral sets of SaS 
laws. If C, is SaS with spectral set Q, then its spectral dual ^* has the 
spectral measure cr*(-) = Sa{Q*,-)- Probabihstic studies of this operation 
call for geometric results concerning pth projection and centroid bodies of 
polar sets. 

Now explore the ordering of SaS vectors based on inclusion relationship 
for their associated star bodies. Write 77 ^ ^ if C for their associated 
star bodies and F^. 

Theorem 10.4. If F^ C F^ for the associated star bodies of SaS random 
vectors ^ and rj with a G (0, 2], then there exist ^ = ^ and fj = rj such that 
\{^,u)\ > \{fi,u)\ a.s. simultaneously for all u G S'^^^ 

Proof. Fix u G S'^^^. Since and {rj^u) are SaS random variables with 

scale parameters \\u\\fi. > 11^11^^, it is possible to define ^ and 77 on the same 
probability space, so that 1(^5^)1 > \{fj,u)\ a.s. By repeating the same ar- 
gument, it is possible to show that finite dimensional distributions of the 
u G S'^"^, are stochastically greater than the finite dimensional dis- 
tributions of \{ri,u)\, u E S'^^^. The statement follows from the continuity of 
the processes, see also [221 Th. 4]. □ 

Definition 10.5. If Fi and F2 are two convex sets representing the unit 
balls in M.'^ with two norms, then the Banach-Mazur distance pBAiiF^i, F2) 
between Fi and F2 (or between the corresponding normed spaces) is defined 
as the infimum of t > such that Fi C AF2 C tFi for an invertible matrix 
A, see [301 Sec. 2.1]. 

The Banach-Mazur distance Pbm{C,v) between two SaS vectors ^ and 
rj with a G [1,2] (needed to ensure the convexity of the associated star 



projections of K onto Hi and H2 respectively. 



□ 
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bodies) is defined as the Banach-Mazur distance between their associated 
star bodies. By Theorem 110.41 Pbm{^,v) is the infimum of t > such that 
C ^ ^ S''^ invertible matrix A. It is well known that the Banach- 

Mazur distance between any d-dimensional space and R"' with the elliptical 
norm is at most \/d. Therefore, for each SaS random vector ^ there exists 
a sub-Gaussian random vector rj with the same characteristic exponent such 
that p_Bj\/(^, rj) < \fd and both r] and ^ can be realised on the same probability 
space as ^ and f] so that 

\{fl,u)\<\{Alu)\<4l\lSi,u)\ a.s. 

holds simultaneously for all u. It is known that (M'^, || ■ ||p) is the farthest 
from the Euchdean among all subspaces of i^p([0, 1]), see [30], Sec. 5.1]. Thus, 
the SaS law with independent components is the farthest one from the sub- 
Gaussian law with the same characteristic exponent. 

Dvoretzky's theorem states that if a natural number n and e > are 
given, then every normed space of sufficiently large dimension d (depending 
on n and e) has an n-dimensional subspace, whose Banach-Mazur distance 
from R'^ with an elliptical norm is less than e. Since section of star bodies 
correspond to projections of SaS vectors, Dvoretzky's theorem imphes that 
each SaS vector with convex associated star body and of sufficiently high 
dimension can be projected onto an n-dimensional subspace, such that its 
projection lies arbitrarily close to a sub-Gaussian law. 

From representation (13. 3p for the characteristic function one immediately 
obtains that if Fi, F2, . . . is a sequence of star bodies corresponding to SaS 
vectors ^1, .^2, • • • with fixed a G (0, 2], then C, (converge in distribution) 
if and only if ^ is an SaS law with associated star body F satisfying ||m||f„ — »• 
\\u\\f as n — s> 00 for all u G M''. It is also possible to provide a version of this 
result for distributions from the domain of attraction of SaS laws. 

Definition 10.6. If is a random vector in M'^, then its associated star 
body at level t is the star-shaped set Ft obtained by (13. 2p using the spectral 
measure at given by 

at(A) = P{-^ e A\ \\ri\\>t} , t>0. 

\m 

The classical limit theorem for convergence to stable random vectors with 
a G (0, 2) implies that if rj belongs to the domain of attraction of SaS law ^ 
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if and only if ||?7|| has a regularly varying tail and at converges weakly to a 
being the spectral measure of ^, see pQ. 

If a > 1, the weak convergence of measures at is equivalent to the Haus- 
dorff convergence of the corresponding L^-zonoids Kt = F^. This can be 
proved in the same way as for a = 1 on [H], p. 184]. For general a G (0, 2], 
the weak convergence of at to a is equivalent to the pointwise convergence 
of norms \\u\\f^ for u G S^^^ together with the convergence of the integrals 
of the norms over the unit sphere. Note that the latter convergence implies 
the convergence of the total masses of at. 

11 James orthogonality 

The associated zonoid can be used as the scale parameter of SaS random 
vector ^ in case a G [1,2]. For general a G (0,2], the star body plays the 
role of the inverse scale parameter. Based on this observation, it is possible 
to generalise several concepts that have been defined only in the univariate 
and bivariate cases or for a > 1. 

The covariation norm frills of SaS random variable rj is defined to be 
the scale parameter of rj, i.e. |||?7|||q, = a if and only if ^n{u) = e~"'^"°, see 
[171 Sec. 2.9]. The family of SaS random variables obtained as linear 
combinations of the coordinates of SaS random vector ^ in M'^ becomes a 
normed space if is equipped with the covariation norm. If ^ has the 
associated star body F and rj = (m,^), then \lri\l^ = \\u\\f, i.e. (5^, ||| ■ |||^) is 
isometric to (M'', || ■ \\f). 

The definition of normality in normed linear spaces goes back to G. Birkhoff 
(1935), see [5(1, Sec. 3.2]. If {X, || ■ ||) is a normed linear space, then x is normal 
to y (notation x H ?/) if ||x + c?/|| > ||x|| for all c G M. This concept was later 
explored by R.C James, and appeared under the name James orthogonality 
in the literature on stable laws. 

If {^1, ^2) are two jointly SaS random variables with a G (1, 2], then ^2 is 
said to be James orthogonal to ,^1 (notation ^2 ^ ^1) if |||c^i + ^2|||q, > 111^2111^ 
for all c G M. The James orthogonality condition can be written as 

i^i^i + «26IIL ^ 1^2! ■ III6IIL ' " = (^1' "2) e . 

If a G (1, 2], we have ^2 H 6 if and only if [^1, ^2]^ = 0, see ^ Prop. 2.9.2]. 
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Theorem 11.1. If (^1,^2) is SaS in with a G (1,2] and the associated 
star body F, then ^2 ^ ^1 if and only if F cM.x [—a, a], where a = Pf((0, 1)). 

Proof. Since the scale parameter of (wi^i + ^2^2) equals h{K,u), the James 
orthogonality condition reads h{K, u) > h{K, (0, ^2)) for all u = {ui, U2) G M. 
By passing to the radial function of F = K* , we see that 

WixW 

If r{6) = pF{cos9,sm9), then 

r{e) < |sin^^|-VF((0,l)), 

which immediately implies the statement, taking into account the equation 
of M X [—a, a] in polar coordinates. □ 

Theorem 111.11 immediately implies that independent SaS variables are 
James orthogonal and that the James orthogonality implies independence in 
the sub-Gaussian case, where F is an ellipsoid. 

The isometry between {S^, ||| ■ |||^) and (R'^, || ■ ||^) makes it possible to 
extend the James orthogonality concept for a G [1,2] (i.e. allow for a = 1) 
and immediately yields the following result. 

Theorem 11.2. Let ^ be SaS with a G [1,2] and associated star body F. 
For each u,v E W^, we have u) H v) if and only if u -\ v in (R"', || ■ ||^). 

Therefore, orthogonality property of SaS random variables from re- 
duces to orthogonality in the normed space (R*^, || ■ H^) if F is convex. It is 
also possible to extend the orthogonality concept for all a G (0, 2] as long as 
F is convex. 

It is known that the orthogonality is symmetric in a normed space of 
dimension at least 3 if and only if the space is Euclidean, i.e. F is an ellipsoid, 
see [Sni Th. 3.4.10]. The corresponding probabilistic result is a part of [ITJ 
Prop. 2.9.3]. In dimension d = 2 the orthogonality is symmetric if and only if 
the boundary of F is a Radon curve, see [B] and [50, p. 94]. The corresponding 
question for SaS laws was posed as an open problem in [47, p. 109]. Recall 
that dF is a Radon curve if and only if the boundary of F in the second 
and fourth quadrants coincides with the boundary of the projection body of 
K = F*. 
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The James orthogonality is a property of the associated star body or 
associated zonoid of an SaS law and is not directly influenced by a. If it 
holds for an SaS law, then it applies for all symmetric stable laws that share 
the same associated star body. 

It is also possible to define multivariate extensions of the James orthog- 
onality concept. 

Definition 11.3. If ^ and 77 are SaS in M'^ with a G [1,2], then 77 is said to 



(i) James orthogonal to ^ (notation r] -\ ^) if the associated zonoid of + 77 

contains the associated zonoid of rj for all c G M; 

(ii) strongly James orthogonal to ^ (notation rj -\s v) if i'^^v) is James 

orthogonal to {u,C,) for all u,v & W^. 

The strong James orthogonality is linear invariant, i.e. all linear trans- 
formations preserve this property. It is easy to see that if rj -\s C,, then 
the associated zonoid of (c^ + rj) contains the associated zonoid of r] for all 
c G M, i.e. the strong orthogonality implies (i). For this it suffices to note 
that this associated zonoid has the support function h{K, {cx,x)) and apply 
Definition lll.3( ii) with u = cx and v = x. 

Theorem 11.4. If {^,r]) is SaS in M?'^ with a G [1,2] and the associated 
star body F, then rj -\s C, if and only if 



for all u = (til, . . . , Mrf, 0, . . . , 0) and v = (0, . . . , 0, Wi, . . . , f^). Further- 
more, T] -\ ^ if and only if llll.l]] holds for u = (c, . . . , c, 0, . . . , 0), v = 
(0, . . . , 0, 1, . . . , 1) and all c G M. 

Proof. For each u',u" G M'^, the scale parameter of c{u',C,) + {u",ri) is 
||(cu', By the condition, this is at least ||(0, which is the scale 

parameter of {u",!]). □ 

If C, and 1] from Theorem 111.41 are independent, then F is the a-star sum 
of the associated star bodies of ^ and 77, i.e. 



be 



(11.1) 



\\u + v\\p = \\u + v\\p^ + \\u + v 
i.e. T] is strong James orthogonal to C,- 
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